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The Zeeman splitting of the donor spectra in cubic- and 
hexagonal- GaN are studied using an effective mass theory ap- 
proach. Soft-core pseudopotentials were used to describe the 
chemical shift of the different substitutional dopants. The 
donor ground states calculated range from 29.5 to 33.7 meV, 
with typically 1 meV higher binding in the hexagonal phase. 
Carbon is found to produce the largest donor binding en- 
ergy. The ionization levels and excited states are in excellent 
agreement with Hall and optical measurements, and suggest 
the presence of residual C in recent experiments. 



Shallow donors in III-N semiconductors are one of 
the key fundamental elements in the fabrication of 
short-wave length optical and high-frequency electronic 
devices. El Despite the great activity and considerable 
progress on the growth quality and characterization of 
these materials, the origin of the n-type conductivity in 
nominally undoped GaN is still not well established. Na- 
tive defects, such as N vacancies, and residual impurities 
like O and Si have been invoked as the prime sources of 
the high electron concentrations (~ 10 17 cm~ 3 ) seen in 
experiments.l3'u However, the donor ionization levels in- 
troduced are not fully known. Temperature dependent 
transport experiments have indicated binding energies of 
about 36 meV for residual donors and 26 meV for Si.u Re- 
cent Hall effect measurements for high-energy electron- 
irradiated GaN/sapphirc samples suggest the appearance 
of a donor state from the N vacancy with a donor energy 
of about 64 meV, much larger than the value expected 
for residual donorsjj Study of electronic excitations by 
Raman scattering of unintentional donors gives energies 
ranging from 18.5 to 30 meV for cubic and hexagonal 
GaN.B Although O is known to display a shallow level 
state, it behaves as a deep donor under high hydrostatic 
pressure (> 20 GPa). However, this is not the case for 
Si.0 

Arrhenius plot analyses of the electron concentration 
of O- and Si- implanted in GaN yields activation ener- 
gies of about 29 meV for both dopants. □ On the other 
hand, magneto-transmission measurements of Zeeman 
splitting of the donor spectra in Si-doped GaNpfilms are a 
bit higher, and consistent with other reports.!] However, 
more recent infrared spectroscopy studies of the Zeeman 
splitting of donors in undoped GaN grown by hydride va- 
por phase epitaxy (HVPE) show Is — 2p± transitions at 
23.3 and 25.95 meV (in the limit of zero magnetic field), 



apparently due to two different donor speciesJj Ground 
states at 31.1 and 33.8 meV, respectively, were inferred 
from these transitions. While the lower transition was 
believed to come from Si, the origin of the larger one was 
unclear. 

Carbon is known to be an amphoteric impurity in GaN, 
but it is widely believed to predominantly introduce ac- 
ceptor levels. Recent quantum-molecular dynamics cal- 
culations by Boguslawski et aZ£3 showed that C substi- 
tuting in the Ga site gives rise to shallow donor states. 
Although the formation energy estimates indicated that 
C atoms do prefer to occupy the N site, there is still the 
possibility that some C atoms will be incorporated into 
the Ga vacancies. In these calculations, C was found 
to have also an excellent solubility, and an efficient self- 
compensatkmof acceptors by carbon donors.E3 Moreover, 
Zhang et a/.til have recently reported photoluminescence 
experiments of C-doped GaN which clearly demonstrate 
that the incorporation of carbon can produce a signifi- 
cant yellow luminescence, and that there is a continued 
C-compensation by the generation of additional shallow 
donors. 

In this work we calculate binding energies for substitu- 
tional donor impurities for different species. Our results 
indicate that residual carbon in GaN films could indeed 
give rise to additional shallow donors with binding en- 
ergy close to those previously assigned to Si and O. To 
our knowledge, the essential effect of central cell correc- 
tions for the different impurities on the donor ionization 
levels in GaN has not been discussed before. We also 
introduce here the effective mass anisotropy in wurtzite 
GaN, which is shown to be important in the calculation 
of impurity states. 

The purpose of this letter is then twofold. First, we 
report effective mass theory (EMT) calculations of the 
binding energy of Si, O and C shallow donors in GaN for 
both crystal phases, cubic (zincblende) and hexagonal 
(wurtzite). In this EMT calculations, the anisotropy of 
the conduction band and the dielectric constant for h- 
GaN are explicitly included. Furthermore, the central 
cell corrections introduced by the different impurities are 
taken into account via soft-core atomic pseudopotentials. 
The second purpose is to examine the Zeeman splitting of 
the excited states of substitutional impurities under the 
presence of an external magnetic field, and to compare 
our estimates with the available experimental results. .— . 

Near the T point (k = 0), the lower conduction bandtZl 
in hexagonal GaN is ellipsoidal, with different longitu- 
dinal and transverse (with the c-axis) effective masses 



1 



(which are clearly equal in the cubic phase). Within ef- 
fective mass theory, and neglecting spin splitting,E3 the 
Hamiltonian for the envelope wave function under an ex- 
ternal magnetic field along the z-axis (parallel to the crys- 
tal c-axis), is given by 
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where we have used a symmetric gauge, fj,j_ is the ef- 
fective Bohr magneton, = ehj2m\c^ and L z is the 
orbital momentum operator along the principal axis (in 
units of h). Here, $ G = hc/e denotes the magnetic flux 
quantum, and U{r) is the impurity-donor potential de- 
fined as the difference between the bare pseudoppten- 
tial of the impurity and that of the host atomJla The 
pseudopetentials we use here are those of Jansen and 
Sankeyta who fitted to.-the functional form of Bachelet- 
Hamann-Schluter fame.E3 These pseudopotentials consist 
essentially of two parts. A local short-range angular- 
momentum-dependent potential Vi(r), and a long-range 
core potential V core (r), where 



V l (r)=J2iMl) + r 2 A l+3 (l)}e- 



,(i)r 2 



and 



V core (r) ^]Taerf[(ar e ) 1/; V] 



(2) 



(3) 



Here, Ai(l), cti(l), Ci and af re are the fitting parame- 
ters to ab initio pseudopotentials, with Z v denoting the 
valence charge. This type of pseudopotential has been 
used widely and with great success in calculations of tha 
electronic and structural properties of semiconductors .cJ 
The extremely good transferability of these pseudopo- 
tentials in the description of atoms in different chemical 
environments has been discussed before by Goedecker et 
aZ.JlJ and is related to the existence of an inert region 
around the nucleus, essentially independent of the chem- 
ical environment. 

The envelope eigenfunctions F(r) of the Hamiltonian 
in Eq. (1) are obtained in a variational approach, and 
chosen here to be an expansion of spherical harmonics 
and a linear combination of hydrogenic radial functions 
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such that F(r) =J2 B ; 
ergy levels are then obtained using the sets {Bj} and 
{/3j} of variational parameters. 

For c-GaN (zincblende), we have used a reliable theo- 
retical value of the isotropic effective mass for the con- 
duction band (m* = 0.19to o )J13 while for the dielectric 
constant we havejised the experimental (and isotropic) 
value of e G = 9.5B [Notice that m* is not known ex- 
perimentally] On the other hand, for the hexagonal 



(wurtzite) phase ot-GaN we have used the experimen- 
tal (and theoretical!!) value of m^j_ p, 0.22 m ,En!j and 
the theoretical value m| = 0.19to o ,E3 for the parallel 
mass. In this case the dielectric constant is obtained from 
the combination e* = (2eJ_7 3 + e7T 1 /3)— 1 = 9.8, where 
e± and arc sct_tp the experimental values of 9.5 and 
10.4, respectively.c3 Due to the moderately polar behav- 
ior of these materials, the effect of the electron-phonon 
interaction on the shift of the energy levels is considered 
through a polaron correction. A Frohlich coupling con- 
stant of qj? = 0.49 has been fixed for both crystal phases. 
This value gives rise to an 8% enhancement of the donor 
states.B 

Figure 1 shows a plot of the position of the electronic 
transitions Is — 2p a and Is — 2p± as a function of mag- 
netic field B for O and C donors in h-GaN. The field is 
parallel to the c-axis direction. The ground state and 
transitions to the first excited states at zero field are 
shown in Table I. For O we find a ground state energy 
of 31.4 meV, and 33.7 for C, with internal transition en- 
ergies Is — 2p± at 23.3 meV and 25.1 meV respectively. 
These values are in excellent agreement to those found 
by Moore et ala at 23.3 and 25.9 meV. Notice we have 
used the same effective mass (m* ± — 0.22 m ) and dielec- 
tric constant (e* = 9.8) as those authors. However, they 
assumed that the lower transition arises from Si donors 
while the higher one was associated with a second donor 
of unknown origin, presumably O. As we show in Table 
I, Si has an EMT ionization energy of 30.4 meV with 
a Is — 2p a transition (which does not shift appreciably 
with field) at 22.8 meV. These results would indicate that 
the binding energy attributed in Ref. ^ to Si is proba- 
bly due to O, and that the second donor may be indeed 
due to residual C-impurities. Since residual C was not 
detected, and the samples were growth by HPVE with 
no carbon sources, it would seem unlikely that carbon 
donors would be present. However, as we mentioned in 
the introduction, small C concentrations of < 2 x 10 16 
cm~ 3 have been found in unintentionally doped samples 
grown by HPVE.O Those concentrations were below the 
detection limit of the secondary ion-mass spectroscopy 
(SIMS), but well within the photoluminescence limits.EJ 
This would indicate that similar SIMS non-detection of 
C impurities in samples obtained by other growth tech- 
niques does not unequivocally prove the absence of small 
quantities of this donor species, allowing our identifica- 
tion to be correct. Carefully monitored C doping to- 
gether with infrared absorption measurements would be 
needed to clarify this point. 

We would like to point out that the transition Is — 2p± 
we find in Si at 22.3 meV, is in good agreement with 



the transition at 21.7 meV observed by Wang et aln by 
infrared absorption spectroscopy of purposely Si-doped 
GaN grown by HVPE. Since the lincwidths observed were 
rather large (~ 5 meV) comparecLwith widths reported 
for nominally undoped samples ,E3 it is likely that the 
small difference of 0.6 meV here is within the uncertainty 
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of the peak position measured. 

Notice also that due to the anisotropy of the effec- 
tive mass and dielectric constant in the wurtzite phase of 
GaN, the 2p a and 2p± states are non-degenerate even at 
zero field. This effect is exhibited in Figure 1 as a shift 
on the transition energy of about 0.5 meV for both im- 
purities. Notice that in Ref. |J a transition at about 17 
meV has also been observed, which remained unaffected 
as the magnetic field was turned on. That feature was 
assigned to the Is — 2p a transition of the donor with 31.1 
meV binding energy, arguing that its origin was related 
to the non-isotropic conduction band of GaN. As we can 
observe in Table I and Figure 1, our EMT calculations 
indicate that such kind of transition does not exist in ci- 
ther donor considered here. Since the associated binding 
energy of such state would be much smaller, the identity 
of that species is not any of the obviously present, and it 
remains a puzzle. 

In Figure 2 we show the transition energy for O and 
C donors as a function of the magnetic field for c-GaN. 
As seen there and in Table I, the ionization levels for this 
phase are close to those of the hexagonal- GaN. We find 
that the difference between the two phases is typically 1 
meV for each donor, being larger in the hexagonal-GaN 
host, partly due the larger effective mass. Notice also 
that the Is — 2p and Is — 2p± transitions are all de- 
generate at zero magnetic field, as one expects in this 
isotropic-band system. Unfortunately, no definitive ex- 
perimental results exist yet for the shallow dnnor levels 
in c-GaN with which to compare our resultso 

Finally, we should mention that in all our EMT cal- 
culations we have employed (were available) the experi- 
mental values of the effective mass and dielectric constant 
for GaN. We find that donor levels are somewhat sensi- 
tive to the variation of these, parameters. For example, 
an uncertainty of ± 0.01m o t3 in the conduction effective 
mass would yield a change of less than ±1.5 meV on the 
energy levels. Notice, however, that all energy levels and 
transitions shown in Table I would shift correspondingly, 
maintaining the relative position for the different impu- 
rities, which helps in the determination of their identity 
(e.g., largest binding energy for C). 

To conclude, let us also mention that we find that Ge 
donors show a spectrum very similar to that found in Si. 
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TABLE I. Ground state Is and lower few internal tran- 
sitions for different impurities in GaN for both crystalline 
phases, zincblende (cubic), and wurtzite (hexagonal). Mag- 
netic field is zero here. All energies are in meV. 
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FIG. 1. Zeeman splitting for O and C shallow donors 
in hexagonal GaN. Magnetic field is applied parallel to the 
c-axis. Anisotropy in the effective mass leads to a shift of the 
Is — 2p transition line at B — 0, while transitions Is — 2p± 
are degenerate. 
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FIG. 2. Zeeman splitting for O and C shallow donors in 
cubic GaN. Magnetic field is applied parallel to a symmetry 
axis. Isotropy of the effective mass in this case yields the same 
value for the Is — 2p and Is — 2p± transitions at B — 0. At 
B / the transitions Is — 2p± split accordingly. 
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